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ABSTRACT 

The  first  important  study  of  fracture  statistics  was  that  of 
Weibull.  His  work  was  based  on  the  tacit  assumption  that  only  the 
component  of  stress  normal  to  a crack  plane  contributes  to  its 
fracture,  and  on  the  use  of  simple  analytical  formulas  for  failure 
probability.  Recent  progress  in  short-term  fracture  includes  the 
use  of  more  refined  fracture  criteria  and  a search  for  better  dis- 
tribution functions  for  the  frequency  of  cracks,  based  on  micro- 
structural  considerations.  Use  of  the  critical  value  of  strain 
energy  release  rate  as  a fracture  criterion  leads  to  improved 
agreement  with  experiment.  Consideration  is  also  given  to  the 
statistics  of  fracture  in  static  fatigue  and  in  dynamic  fracture. 


INTRODUCTION' 

Ceramic  materials  are  generally  brittle,  and  nominally  identi- 
cal specimens  of  a brittle  material  may  exhibit  large  variations  in 
fracture  stress,  especially  if  the  specimens  are  small.  When  brit- 
tle materials  are  employed  in  practical  structures,  the  designer 
must  be  able  to  assure  himself  of  an  acceptably  low  probability  of 
failure  during  service.  Typically,  on  the  basis  of  laboratory  data 
on  a limited  number  of  specimens  uniformly  loaded  in  simple  tension 
or  pure  bending,  the  probability  of  failure  must  be  calculated  for 
structural  members  of  different  sizes  and  shapes  and  under  com- 
pletely different  loading  conditions.  The  tool  for  accomplishing 
this  is  fracture  statistics. 


The  first,  important  contribution  to  this  subject  was  made  by 
Keibull  in  1939  [1,2],  and  his  theories  are  still  the  basis  for 
most  calculations  in  this  field.  He  exploited  the  analogy  between 
a stressed  brittle  structure  and  a loaded  chain,  which  breaks  when 
the  strength  of  its  weakest  link  is  exceeded.  Before  discussing 
weakest  link  theory  [WLT],  it  is  appropriate  to  enquire  into  its 
range  of  applicability. 

Clearly,  the  fracture  of  a single  fiber  is  similar  to  the 
breaking  of  a chain.  The  fiber  may  be  regarded  as  a linear  array 
of  very  short  elements,  and  the  fracture  o*  any  element  causes 
failure  of  the  fiber  as  a whole.  Similarly,  when  an  isolated 
crack  in  an  elastic  body  is  loaded  normal  to  its  plane,  it  will 
become  unstable  and  grow  catastrophically,  causing  fracture.  The 
fracture,  stress  of  the  entire  body  is  that  of  the  weakest  crack. 

On  the  other  hand,  there  are  many  situations  to  which  WLT  does  not 
apply.  For  instance,  in  a bundle  of  fibers  the  first  fiber  fail- 
ure does  not  ordinarily  cause  failure  in  the  bundle  as  a whole; 
rather,  failure  is  the  result  of  damage  accumulation.  Similarly, 
in  some  stress  states  (e.g.  pure  compression)  a crack  in  an  elastic 
body  will  generally  grow  in  a direction  which  results  in  crack 
arrest.  Here,  also,  failure  is  a result  of  progressive  damage, 
and  WLT  does  not  apply. 

Ideally,  fracture  statistics  should  be  based  on  a proper  con- 
sideration of  three  elements  — extreme  value  statistics,  fracture 
mechanics,  and  material  microstructure.  Weibull's  theory  is  based 
almost  exclusively  on  the  first  element.  In  recent  years,  progress 
has  been  made,  in  incorporating  the  other  two  elements  into  weakest 
link  theory.  Tne  purpose  of  the  present  paper  is  to  outline  the 
current  status  of  fracture  statistics  and  indicate  some  directions 
of  expected  future  progress.  In  the  effort  to  give  a clear  and 
concise  description  of  the  field,  historical  perspective  has  been 
somewhat  slighted.  The  primary  emphasis  is  on  the  underlying 
physical  concepts  rather  than  computational  techniques.  It  is 
hoped  that  those  whose  contributions  may  have  been  omitted  or 
underemphasized  will  understand  and  forgive. 


WEAKEST  LINK  THEORY  (WLT) 

We  first  derive  the  fundamental  equation  of  weakest  link 
theory.  Let  it  be  assumed  that  a stressed  solid  can  fail  due  to 
any  of  a number  of  independent  and  mutually  exclusive  mechanisms 
or  causes,  each  involving  infinitesimal  probability  of  failure 
(APf)j^.  The  probability  th.it  the  i'th  mechanism  will  not  cause 
failure  is  (Ps) i = 1 - (APf)f.  The  overall  probability  of  sur- 
vival is  the  product  of  the  individual  probabilities  of  survival 
i.e . , 
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ps*  II  (ps).  • IX  {l  - ("'f)  } 

iii  i; 

s TIoxp  - (AVf)J=exp  F-  £ (APf).  • (D 

Tl;e  sum  of  the  individual  probabilities  of  failure  appearing  in  the 
final  equality  above  was  called  by  Weibull  the  "risk  of  rupture" 
and  was  given  the  symbol  B. 

In  evaluating  B,  Weibull  added  together  the  probabilities  of 
failure  of  all  the  elements  of  volume  AV  in  the  entire  body.  The 
probability  of  failure  of  the  i ' th  element  AV.  in  simple  tension 
a,  for  example,  is  1 

(APf)  = n(o)  AV.  , (2) 


where  n(a)  is  the  number  of  flaws  per  unit  volume  with  a strength 
less  than  a.  If  n(o)  is  less  than  unity,  it  can  be  regarded  as  the 
probability  that  such  a flaw  will  occur  in  a unit  volume.  Thus, 
the  probability  of  failure  is  given  by 

Pf  = 1 " PS  = 1 “ CXp  ["  / dV  n(0)]  * (3) 

For  a uniformly  stressed  body 

Pf(a)  = 1 - exp  [-  V n(a)]  , (4) 

or 

n(a)  - “ s,n  Pg-1  . (5) 

If  Pf (a)  is  determined  by  testing  a number  of  specimens  in  simple 
tension,  F.q.  (5)  can  be  used  to  obtain  n(o)  and  Eq.  (3)  can  then  be 
used  to  evaluate  the  probability  of  failure  of  a body  of  arbitrary 
size  and  shape  and  nonuniformly  distributed  tensile  stress.  By 
assuming,  as  Weibull  did,  that  compressive  stresses  do  not  contrib- 
ute to  fracture,  the  above  procedure  can  also  be  used  to  calculate 
the  probability  of  failure  in  pure  bending,  or  indeed  any  combina- 
tion of  uniaxial  tension  and  compression.  A limitation  to  this 
procedure  is  that  in  ignoring  fracture  due  to  compression,  it  does 
not  adequately  account  for  failure  in  all  stress  states. 
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Weibull  also  gave  a procedure  for  calculating  Pf  for  polyaxial 
stress  states;  when  the  failure  statistics  for  simple  tension  are 
known.  Basically,  his  procedure  is  to  calculate  B by  averaging 
the.  tensile  stress  in  all  directions  [1,3].  This  is  intuitively 
plausible,  but  not  rigorous,  and  it  has  been  questioned  by  some 
investigators.  Barnett  et  al.  [4]  have  formulated  a number  of 
alternative  procedures  they  considered  equally  plausible.  Presum- 
ably, in  part  at  leant,  because  of  doubts  concerning  Weibull's 
procedure,  some  investigators  [5,6,7]  have  chosen  to  use  an 
approximation  in  which  it  is  assumed  that  with  respect  to  fracture, 
the  principal  stresses  act  independently,  i.e., 

W °2'  °3>  ■ V°l>  W • (6) 


where  0^ , 0^,  are  the  principal  stresses. 

I 

There  is  some  evidence  that  Weibull  had  second  thoughts  con- 
cerning the  validity  of  his  recommended  procedure,  since  he  wrote 
in  1966  [8]:  "Another  problem  of  a more  theoretical  nature  will  be 

to  deduce  the  effect  of  bi-  and  triaxiality  on  the  distribution 
functions  of  one- dimensional  stresses.  If  the  principal  stresses 
are  acting  independently  of  each  other,  and  it  seems  that  such 
materials  may  exist,  then  the  problem  may  be  soluble  along  lines 
previously  sketched.  In  other  cases  the  solution  is  very  intricate 
and  will  certainly  require  close  examination  of  the  physical 
behavior  of  the  material  in  question." 

WLT  FOR  POLYAXIAL  STRKSS  CONDITIONS 

To  put  weakest  link  theory  for  polyaxial  stress  states  on  a 
firm  physical  foundation,  we  will  make  the  explicit  assumption  that 
the  flaws  responsible  for  fracture  are  microcracks  in  the  material. 

We  will  further  assume  that  the  cracks  do  not  interact,  and  that 
each  crack  has  a critical  stress  0C  defined  as  the  remote  tensile 
stress  applied  normal  to  the  crack  plane  which  will  cause  fracture. 

Fracture  under  combined  stresses  occurs  when  the  effective  stress 
0C  acting  on  the  crack  is  equal  to  Oc.  The  effective  stress  is 
some  function  of  the  applied  stresses  at  the  location  of  the  crack, 
and  its  precise  form  depends  on  the  fracture  criterion  employed. 

The  search  for  fracture  criteria  leads  us  to  fracture  mechanics. 

Before  discussing  particular  fracture  criteria,  however,  we  formu- 
late the  theory  of  fracture  statistics  for  an  arbitrary  fracture 
criterion  by  leaving  unspecified. 

The  potential,  causes 
For  purposes  of  analysis. 
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of  failure  are  the  individual  cracks, 
it  is  convenient  to  group  the  cracks 
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according  to  location,  the  applied  stress  state,  and  crack  critical 
stress.  We  assume  that  the  stress  state  varies  slowly  so  that 
within  a volume  element  AV  all  cracks  will  be  subject  to  the  same 
macroscopic  stress.  We  also  assume  that  the  material  is  uiacro- 
scopically  homogeneous,  so  that  a function  N(cc)  can  be  defined  as 
the  number  of  cracks  per  unit  volume  having  a critical  stress  equal 
to  or  less  than  oc . The  probability  that  a crack  having  a critical 
stress  in  the  range  Oc  to  ac  + doc  exists  in  volume  element  AV  is 


the  AV  [dN(a  )/dO  1 da  . 

c c ‘ c 


If  such  a crack  is  actually  present,  the  probability  that  it 
will  fracture  depends  on  its  orientation,  the  stress  state,  and 
the  fracture  criterion.  We  assume  there  is  a solid  angle  ft  such 
that  fracture  of  a crack  will  occur  if  and  only  if  its  normal  lies 
within  ft.  This  means  that  if  the  normal  lies  within  ft,  oe  > 0C  , 
where  Oe  is  the  effective  stress  corresponding  to  the  fracture  cri- 
terion selected.  If  the  cracks  are  randomly  oriented,  the  proba- 
bility that  a crack  will  fracture  under  the  applx€:d  stress  Z is 

ft(Z,0  ) / 4 if . 
c 

Now  the  probability  of  failure  due  to  a crack  in  the  critical 
stress  range  rtoc  located  in  volume  element  AV  if;  the  product  of  the 
.above  probabilities,  i.e.. 


/ dN(o  ) 


'fi(Z,a  ) 

/nr 


(/) 


Substituting  Eq.  ('/)  into  Eq.  (1),  and  changing  sums  into  integrals, 
we  obtain 
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?S  = exp  [-  / dV  / dac  f-]  . (8) 

We.  note  in  passing  that  N(ac)  is  independent  of  stress  state 
and  depends  only  on  the  material . Also,  since  oc  is  defined  ns 
the  stress  which  causes  fracture  when  applied  normal  to  the  crack 
plane,  N(cc)  can  be  converted  directly  into  crack  size  distribution 
when  Kjc  is  known.  In  these  respects  it  differs  from  Weihull’s  n 
and  the  g(S)  of  McClintock  and  others  [9,10,11  ].  These  functions 
represent  the  number  of  cracks  per  unit  volume  that  will  be  frac- 
tured by  a particular  applied  stress  state,  and  derivative  of  such 
a function  respectively,  and  they  depend  on  the  stress  ratio.  For 
the  particular  case  of  hydrostatic  tension,  N(c0  - n (cl) and  K1  (a)  = 
g(a) , because  then  ft  - 4tt  or  zero  depending  on  the  stress  level. 


In  simple  tension  and  equibiaxial  tension,  analytical  expres- 
sions can  be  found  for  ft  for  some  fracture  criteria,  at  least,  and 
used  together  with  Eq.  (3)  to  evaluate  l*s.  In  the  general  case, 
we  find  ft  by  integrating  dft  over  the  range  in  which  oc  > 0C.  One 
way  of  accomplishing  this  is  to  integrate  over  the  entire  angular 
range  but  include  a suitable  operator  li  in  th.  integral: 


PS  " 


exp 


in 


dv  dft  da  II (a  ,o  ) 

c e c da 

c 


(9) 


where 


H(a  ,o  ) = 1 when  a > a 
e c e c 

= 0 where  a < a . (10) 

e c 

ow  carry  out  the  integral  over  ac  first,  with  the  result 

Pg  = exp  \ - jj  dV  dft  N(ae)]  . (11) 

The  effective  stress  causing  fracture  is  a function  of  both 
the  component  of  stress  normal  to  the  crack  plane,  an,  and  the 
shear  stress  T parallel  to  the  crack  plane.  An  approximation  fre- 
quently employed  is  the  assumption  that  the  cracks  are  shear- 
insensitive,  i.e.,  that  ae  - 0n . Using  this  approximation, 

Eq.  (11)  becomes  equivalent  to  Keibull's  rule  for  polyaxi.nl  stress 
states  [ 1 ] , and  Eq.  (8)  reduces  to  the  equation  of  Batdorf  and 
Crose  [l2],  Eq.  (8)  and  Eq . (11)  are  equivalent  formulations  of 
the  same  theory.  The  physical  justification  of  the  former  is  the 
more  readily  apparent,  but  the  latter  is  more  convenient  for  com- 
putational purposes. 


RELATION  TO  GRIFFITH  THEORY 

In  1924,  Griffith  f 1 3 ] published  a theory  for  the  fracture  of 
solids  under  biaxial  stress  conditions.  He  assumed  the  presence  of 
a largo  number  of  identical  cavities  in  the  form  of  nearly  flat 
elliptical  cylinders  (Griffith  cracks)  with  their  axes  normal  to 
the  plane  of  the  stresses,  but  otherwise  randomly  oriented.  He 
calculated  the  maximum  peripheral  tensile  stress  on  the  free  sur- 
face of  such  a crack  as  a function  of  crack  orientation,  and 
assumed  fracture  would  occur  whenever  this  maximum  value  exceeded 
the  intrinsic  strength  of  the  material. 


Tin's  fracture  criterion  led  to  the  failure  envelopes  shown  in 
Fig.  1.  In  the  tension-tension  quadrant  and  most  of  the  tension- 
compression  quadrant , failure  obeys  a maximum  tensile  stress  cri- 
terion. Fracture  initiates  in  a crack  whose  plane  is  normal  to  the 
largest  principal  tensile  stress.  Thus,  only  the  normal  stress  On 
acts  on  the  crack  plane,  and  tensile-type  fractures  result.  In  the 
portion  of  the  tension-compression  quadrant  near  the  compression 
axis,  material  strength  is  locally  exceeded  when  the  crack  normal 
is  inclined  at  somewhat  more  than  45°  to  the  axis  of  maximum 


a 


2 


vYLT,  shear- 
insensitive 
cracks 


Y/LT,  shear- 

sensitive 

cracks 


Fig.  1.  Comparison  of  Fracture  Theories  for  Biaxial  Stresses 
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tension.  Here  both  0n  and  T contribute  to  failure  and  the  result- 
ing oblique  fracture  is  often  called  a shear-type  failure. 

In  any  crack-based  statistical  theory  it  is  most  unlikely  that 
the  weakest  crack  will  be  oriented  at  exactly  the  angle  of  maximum 
vulnerability,  so  Griffith's  curves  must  be  regarded  as  the  lower 
limit  for  fracture,  i.e.,  they  represent  the  coi.tour  for 
Pf(01,02)  <<:<  1.  A typical  Pf(oi,02)  = 0.5  contour,  for  the  case 
of  shear- insensitive  cracks,  is  shown  as  a dashed  curve.  This 
curve  covers  both  quadrants  almost  completely,  but  does  not  provide 
for  failure  due  to  pure  compression.  In  the  case  of  shear-sensitive 
cracks,  the  entire  Pf(0l,O2)  = 0.5  contour  predicted  by  WLT  will 
fall  slightly  outside  of  Griffith's  failure  envelope  as  shown  sche- 
matically in  Fig.  1.  For  high  compressive  stresses,  the  results 
differ  significantly  from  those  for  shear-insensitive  cracks. 
However,  as  noted  earlier,  in  this  region  WLT  does  not  apply. 

Thus,  Weibull  theory  and  the  more  refined  theories  to  be  discussed 
later  are  limited  to  cases  in  which  no  principal  compressive  stress 
exceeds  the  maximum  principal  tensile  stress  in  absolute  value  by 
more  than  a factor  of  about  3.  This  limitation  is  probably  not  too 
serious  from  a practical  point  of  view  because  the  percent  disper- 
sion in  fracture  stress  is  generally  much  less  in  compression  than 
in  tension  [14,15]. 


CRACK  DENSITY  FUNCTIONS 

The  derivative  of  the  crack  density  function,  N'(0C),  gives 
the  number  of  cracks  per  unit  volume  per  unit  range  in  critical 
stress.  One  might  at  first  be  inclined  to  expect  this  function  to 
be  Gaussian  or  near-Gaussian . In  actuality  it  turns  out  that 
Gaussian  distributions  are  incompatible  with  WLT.  This  is  readily 
demonstrated  as  follows. 

We  have  seen  that  any  WLT  must  take  the  form 

Pg(o)  = - exp  [-  f(o)]  , (12) 

whence 

Pf(0)  = " Ps(0)  = f’(o)  exp  f(a)l  • (13) 

If  the  distribution  is  Gaussian, 

P j(o)  = A exp  [-a  (0  - b) ^ 3 


(14) 


9 


Eq.  (13)  and  Eq.  (14)  are  incompatible,  because  if  f’(0)  is  a con- 
stant, f (o)  cannot  be  a quadratic  function.  Another  reason  for 
excluding  the  Gaussian  distribution  is  that  it  implies  a nonvanish- 
ing probability  of  tensile  failure  when  the  applied  stress  is 
compression . 

Weibull  introduced  two  distribution  functions 

Pf (0)  = 1 - exp  [-  V k0m]  (15) 

and 

P (0)  = 1 - exp  [-  Vk  (a  - o )m]  (o  > a ) (16a) 

f L u u 

=0  (0  < 0u)  (16b) 

Each  may  be  regarded  as  a skewed  Gaussian  distribution  with  a 

skewness,  which  can  be  either  positive  or  negative,  determined  by 

the  value  of  m.  The  2-parameter  form  allows  failure  to  occur  at 

any  positive  value  of  the  tensile  stress,  while  the  3-parameter 

form  implies  that  fracture  cannot  occur  for  a < o . 

u 

The  function  Pf(o)  is  conventionally  determined  by  conducting 
N tests  and  numbering  the  observed  fracture  stresses  CTi...0fj  in 
ascending  order.  It  is  then  usually  assumed  that 


Pf(cV  N + 1 


A more  sophisticated  statistical  treatment  [16]  leads  to  the  con- 
clusion that 


Pf(cV  N + o!a 


The  difference  between  Eq.  (17)  and  Eq . (18)  becomes  insignificant 
for  large  values  of  N and  for  simplicity  we  shall  use  Eq.  (17). 

A simple  technique  for  determining  the  parameters  m,  o , and 
k is  to  write  Eq.  (16a)  in  the  form  u 

£n  £n  (1  - P^)  ' = £n  Vk  + m £n  (a  - Oy)  . (19) 

Next,  £n  £n  (1  - Pf)  1 is  plotted  against  £n  (o  - Ou)  for  various 
assumed  values  of  0^.  The  value  of  adopted  is  that  for  which 


the  N data  points  most  nearly  lie  on  a straight  line.  The  slope  of 
the  line  chosen  is  m,  and  Vk  is  the  value  of  £n  (1  - Pf)~l  for 
0 = 0U  + 1.  When  Weibull's  2-parameter  form  is  used,  0 is  arbi- 
trarily taken  to  be  zero.  u 

Before  passing  on,  we  note  that  a least  squares  fit  in 
in  £n  (1  - Pf)“l  vs.  £n  (a  - au)  space  is  in  general  not  a least 
squares  fit  in  Pf  vs.  O space.  Consequently,  the  Weibull  form 
found  in  the  manner  described  is  not  the  best  possible  fit  to  the 
data.  However,  if  the  scatter  in  the  data  is  not  large,  it  will 
be  adequate  for  most  practical  purposes. 

Fracture  statistics  are  relatively  simple  when  Weibull's 
2-parameter  form  is  used.  If  the  probability  of  failure  in  simple 
tension  is 


Pf  = 1 - exp  (-  V kT  am)  , (20) 

then  that  in  bending  can  be  shown  to  be  [1] 

Pf  - 1 ~ exp  [-  V kfi  0m)]  , (21) 

where  kjj/kx  depends  on  the  cross  section.  For  a rectangular  cross 
section 


kR  = kT/(2m  + 2)  . (22) 

Under  the  uniform  principal  stresses  0^,  0^,  O-,  the  probability  of 
failure  is  given  by 


The  constant  k (02/01,  03/^)  is  an  analytical  function  of  the 
stress  ratios  [17].  Many  other  cases  are  treated  in  [18]. 

Weibull's  3-parameter  representation  has  a similar  universal 
applicability  when  attention  is  limited  to  uniaxial  stress  prob- 
lems. In  going  from  uniaxial  to  polyaxial  stress  states,  however, 
the  functional  form  changes  [19],  i.e.  it  is  not  possible  to  express 
both  uniaxial  and  polyaxial  stress  data  in  the  form 
Pf  = 1 - exp  [-  V k (0  - ou)m].  To  assist  those  desiring  to  use 
Weibull's  3-parameter  representation  in  conjunction  with  his  theory 
for  polyaxial  stress  states,  Dukes  [20]  has  carried  out  parametric 
calculations  with  the  aid  of  a computer. 


Weakest  link  theory  implies  a volume  effect  that  is  sometimes 
misapplied.  Consider,  for  instance,  the  situation  illustrated 
schematically  in  Fig.  2.  Let  us  assume  specimens  were  tested,  so  9 
we  know  Pf(o)  over  the  range  0.1  < Pf  < 0.9.  The  known  region  is 
shown  as  a solid  curve,  and  the  dotted  extension  is  an  extrapola- 
tion of  the  data,  accomplished  by  assuming  that  the  mathematical 
function  used  to  represent  the  test  data  continues  to  be  valid  out- 
side the  range  in  which  it  was  tested.  All  that  WLT  can  tell  us  is 
that  a specimen  10  times  larger  will  have  1C  times  the  risk  of  rup- 
ture over  the  stress  range  tested.  It  does  not  answer  the  generally 
more  interesting  question,  what  stresses  correspond  to 


WLT 


DATA 


SPECULATION 


Fig.  2.  Data,  WLT  inference 

concerning  volume  effect,  and 
speculative  extrapolations. 


probabilities  of  failure  of  the  larger  specimen  in  the  probability 
range  tested,  i.e.,  0.1  < Pf  < 0.9,  unless  the  validity  of  the 
extrapolation  to  lower  stress  levels  can  be  established. 

This  raises  the  question  whether  the  Weibull  distribution 
functions  have  a fundamental  significance  that  might  justify  confi- 
dence in  such  extrapolation.  Weibull  made  no  such  claim;  he 
regarded  them  only  as  convenient  mathematical  tools  having  wide 
applicability  [ 2 1 j . It  is  true  that  they  coincide  with  the  third 
asymptotic  form  of  extreme  value  theory.  However,  the  founders  of 
this  theory  did  not  assert  that  all  distributions  must  approach  a 
limiting  form,  but  only  that  if  a distribution  does  approach  a 
limiting  form,  it  must  be  one  of  the  three  identified  in  the  theory 
of  extreme  values  [22,23] . Basically,  the  limiting  form  results 
when  a distribution  can  be  represented  as  a power  series.  When 
this  happens,  only  the  leading  term  need  be  considered  for  large 
samples  (i.e.  large  specimens  with  very  many  cracks).  We  shall 
show  later,  however,  that  microstructural  considerations  suggest 
that  such  a power  series  is  not  possible,  and  that  as  a result  no 
limiting  form  exists  for  fracture. 


FRACTURE  CRITERIA 

The  direct  stress  applied  normal  to  a crack  plane,  on,  results 

in  a very  high  local  stress  at  the  root  of  the  crack,  while  direct 

stresses  in  the  plane  of  the  crack  do  not.  The  shear  stress  T 

applied  parallel  to  the  crack  plane  also  results  in  very  high  local 

stresses.  Consequently,  the  effective  stress  C causing  fracture 

is  a function  of  both  O and  T. 

n 

The  statistical  analysis  of  fracture  is  appreciably  simplified, 
however,  by  the  assumption  that  Oe  = On,  i.e.  that  the  cracks  are 
shear-insensitive.  In  this  case  the  properties  of  the  cracks  are 
completely  characterized  by  0C,  the  remote  stress  normal  to  the 
crack  plane  that  will  cause  fracture.  There  is  then  no  need  to 
specify  the  size  or  shape  of  the  crack  or  Poisson's  ratio.  And 
because  any  given  crack  is  weakest  when  oriented  normal  to  the 
largest  principal  tensile  stress,  as  pointed  out  in  the  section 
"Relation  to  Griffith  Theory,"  the  approximation  Cfe  = 0r.  is  not  a 
bad  one.  This  approximation  was  explicitly  made  by  Batdorf  and 
Crose  [24],  and  is  implicit  in  Weibull's  treatment  of  polyaxial 
stress  states  [1]. 

The  fracture  criterion  used  by  Griffith  to  take  into  account 
the  effects  of  shear  is  the  assumption  that  fracture  occurs  when 
the  local  tensile  stress  at  some  point  on  the  crack  surface 
exceeds  the  intrinsic  strength  of  the  material.  Many  authors  have 
discussed  the  stress  distribution  around  cavities  of  various  types 


under  various  loading  conditions.  For  present  purposes,  the  most 
concise  and  general  treatment  is  probably  that  of  Mirandy  and  Paul 
[25]  who  give  an  explicit  formula  for  the  maximum  stress  on  the 
surface  of  an  ellipsoid  with  semiaxes  a > b » c under  arbitrary 
applied  loads.  They  find  that 

a ="  + \la2  + t2  F ! , (24) 

max  cE/nvn  J 


where  E is  an  elliptic  integral  that  depends  on  b/a  and  F is  a com- 
plicated function  of  the  geometry  of  the  crack  and  its  orientation 
relative  to  the  shear  stress.  We  have  defined  the  effective  stress 
to  be  in  the  absence  of  shear,  and  some  function  of  an  and  T 

when  shear  is  present.  Accordingly,  we  conclude  that  for  a mate- 
rial containing  cracks  of  only  a single  plan  form  (i.e.  a fixed 
value  of  b/a) 

- I { + Van  + T2  F j ' <25> 

Using  the  formulas  in  [25],  it  can  be  shown  that  for  a Griffith 
crack  (G.C.)  with  its  axis  normal  to  the  applied  shear  stress 


(G.C.) 


while  for  a penny-shaped  crack  (P.S.C.) 
I 


0 

e 


+ t2/(1  - 0.5v)2 


(P.S.C.) 


(26) 


(27) 


where  v is  Poisson's  ratio.  The  latter,  possibly  less  familiar 
result,  is  also  given  in  [26]. 


There  are,  however,  grounds  for  doubting  that  maximum  local 
tensile  stress  represents  an  acceptable  fracture  criterion.  Con- 
sider, for  instance,  an  ellipsoidal  cavity  with  principal  axes 
a ■ 4b  loaded  in  tension  parallel  to  the  c axis.  The  maximum  local 
stress  occurs  in  the  equatorial  plane,  and  it  is  uniform  around  the 
entire  circumference.  Thus,  according  to  the  maximum  local  tensile 
stress  criterion,  all  points  on  the  equatorial  belt  are  equally 
likely  to  fracture  and  the  crack  should  advance  in  all  directions 
in  its  own  plane.  On  the  other  hand,  the  region  which  is  above 
some  fixed  fraction  of  the  maximum  stress  is  four  times  as  wide  at 
the  ends  of  the  b axis  as  at  the  ends  of  the  a axis.  Thus,  if 
stressed  volume  is  important,  growth  should  start  first  at  the  ends 


of  the  b axis.  Similarly,  the  stress  intensity  factor  is  largest 
at  the  ends  of  the  b axis,  where  it  is  twice  as  large  as  at  the 
ends  of  the  a axis.  If  stress  intensity  governs  fracture,  the  b 
axis  should  extend  first,  and  the  crack  should  circularize.  It  is 
generally  agreed  that  this  is  what  actually  happens. 

Another  fracture  criterion  is  the  assumption  chat  fracture 
occurs  when  the  elastic  strain  energy  released  when  the  crack  grows 
is  equal  to  the  energy  stored  in  the  newly  created  free  surfaces, 
i.e.  when  the  strain  energy  release  rate  reaches  some  critical 
value.  The  strain  energy  release  rates  are  well-established  for  a 
Griffith  crack  or  penny-shaped  crack  growing  in  its  own  plane. 

For  the  former  [27] 


c ^ [K?  + Kn]  • (28) 

where 

E = elastic  modulus 

Kx  --  a /0a~  (29) 

. K = T /iwT  . (30) 

Thus, 

a =\/a2  + t2  (G.C . ) . (31) 

e V n 

For  the  latter  [28] 


G = (1  V2) 

[k2t  + k2+(  1 \ 

K2  1 , 

(32) 

E 

where 

L 1 11  \i  _ v2/ 

III  ’ 

Kt  = 2o 
I n 

s/a/iT 

(33) 

*n  ' 1 T 

v/a/ir  sin  y/ (2  - v) 

(34) 

Km  ■ 4 T 

Ja/v  cos  y (1  - v)/(2 

- v) 

(35) 

The  minimum  value  of  G 

occurs  at  the  points  on 

the  crack  periphery 

where  only  modes  I and  II  are  involved,  and  here 
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a *=  \/o 2 + T2/(l  - 0 .5v) 2 (P.S.C.)  . (36) 

e V n 

Unfortunately,  under  the  combined  action  of  on  and  T,  cracks 
do  not  extend  in  their  own  plane  so  the  above  expressions  are  not 
strictly  correct.  There  is  at  present  no  consensus  regarding  the 
proper  fracture  criterion  for  cracks  extending  out  of  their  plane. 
Swedlow  [29]  listed  nearly  thirty  differing  treatments  of  this  sub- 
ject that  had  appeared  by  1975,  and  others  have  been  proposed  since 
that  time.  We  will  not  attempt  to  employ  any  of  these  theories  in 
this  review.  Instead,  we  will  use  the  equations  just  derived  not- 
ing, however,  that  the  true  value  of  oe  must  be  smaller  than  that 
found  for  coplanar  extension.  That  is  because  the  lowest  instabil- 
ity mode  is  the  mode  that  actually  occurs. 

Some  idea  of  the  differences  between  the  various  fracture  cri- 
teria discussed  herein  can  be  obtained  by  assuming  the  material 
obeys  the  Weibull  2-parameter  form  [Eq.  (15)],  and  plotting  the 
ratio  of  the  risk  of  rupture  under  equibiaxial  tension  to  that  for 
uniaxial  tension  vs.  the  parameter  m.  This  is  done  in  Fig.  3 
which  is  adapted  from  [17]  and  [30].  The  ratio  is  largest  for 
shear  insensitive  cracks.  The  remaining  curves  are  for  the  cri- 
'teria  in  Eqs.  (26),  (27),  (31),  and  (36)  respectively.  We  note 
that  the  ratio  kQC]  hi as^uniax  increases  with  m for  all  criteria, 
and  is  smallest  for  the  fracture  criteria  based  on  strain  energy 
release  rate.  Also  included  in  the  comparison  is  the  assumption 
that  with  respect  to  fracture,  the  principal  stresses  behave  inde- 
pendently. This  leads  to  a constant  ratio  of  2.  The  comparison 
suggests  that  for  m > 3 the  independence  assumption  will  lead  to 
unconservative  estimates  of  the  statistics  of  failure  in  biaxial 
tension.  We  note  in  passing  that  for  the  low  values  of  Pf  desired 
in  most  structural  applications,  the  failure  probability  ratio  is 
the  same  as  the  k ratio. 

Experimental  data  on  uniaxial  and  equibiaxial  bending  of  alu- 
mina plates  were  recently  obtained  by  Sines  and  Giovan  [30],  In 
Fig.  4 the  equibiaxial  data  are  approximated  by  a Weibull 
2-parameter  curve.  The  uniaxial  results  were  computed  for  shear- 
insensitive  cracks  and  also  for  the  fracture  criteria  listed  in 
Eq.  (27)  and  Eq.  (36).  As  would  be  anticipated  on  the  basis  of  the 
preceding  discussion,  the  agreement  is  best  using  Eq.  (36).  For 
details  of  the  analysis,  see  [30]. 

The  relative  merits  of  various  fracture  criteria  can  also  be 
tested  in  other  ways.  For  instance,  Petrovic  and  Mendiratta  [32] 
tested  the  variation  of  tensile  fracture  stress  of  controlled  sur- 
face cracks  with  crack  angle.  The  conclusion  from  this  experiment 
is  that  Eqs.  (31)  and  (36)  are  in  reasonably  good  agreement  with 
experiment,  and  are  the  best  of  the  criteria  listed  in  this  paper. 
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Fig.  3.  Relation  Between  Failure  Probability  Under  Equibiaxial 
Loading  and  That  for  Uniaxial  Loading  as  a Function  of 
Weibull  Parameter  m. 


WLT  FOR  ANISOTROPIC  MATERIALS 

Anisotropic  brittle  materials  may  be  anisotropic  with  respect 
to  their  elastic  properties,  their  fracture  properties,  or  both. 

The  elastic  anisotropy  will  influence  the  stresses  resulting  from 
a given  load  application,  but  will  not  directly  affect  the  calcu- 
lation of  probability  of  failure,  which  depends  only  on  the  stress 
state  and  crack  distribution. 

Two  different  approximate  techniques  have  been  proposed  for 
determining  the  fracture  statistics  of  certain  types  of  polygraph- 
ite that  are  isotropic  in  one  plane  but  have  a lower  fracture 
stress  in  the  direction  normal  to  the  plane  of  isotropy.  In  one 
approach,  the  cracks  were  assumed  randomly  oriented,  but  they  were 
given  critical  stresses  that  varied  with  orientation  [ 3 1 3 - In  the 
other,  the  critical  stresses  were  assumed  unaffected  by  orientation, 
but  cracks  were  given  a preferred  orientation  [32].  It  is  likely 
that  there  is  actually  both  a preferred  orientation  and  a variation 
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Fig.  4.  Theory  Compared  with  Data  of  Giovan  and  Sines 


in  crack  strength  with  orientation,  but  a theory  incorporating  both 
would  be  considerably  more  complex  ana  not  really  needed,  at  least 
for  graphite.  Both  theories  that  attribute  the  entire  anisotropy 
to  only  one  effect  appear  to  be  in  satisfactory  agreement  with  the 
available  data. 


WLT  FOR  SURFACE  DISTRIBUTED  CRACKS 

In  his  first  paper  on  the  statistics  of  fracture,  Weibull 
included  a discussion  of  surface-distributed  flaws  [l].  In  this 
treatment,  the  only  change  was  to  replace  B = J°  n^(a)  dV  by 
B *=  f n 2 (a)  dA  where  A represents  area,  and  n2  might  or  might  not 
be  the  same  as  nj.  Thus  the  analogy  to  the  theory  for  volume  dis- 
tribution flaws  is  very  close. 

In  a more  recent  analysis  the  fracture  statistics  of  surface 
distributed  cracks  have  been  worked  out  on  the  assumption  that  crack 
planes  are  always  normal  to  the  material  surface  [35].  As  a result, 
the  orientation  of  a crack  is  specified  by  a single  angle,  whereas 
for  volume  distributed  cracks  two  angles  are  required.  The 
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analogue  to  Eq.  (8)  applicable  to  surface  distributed  cracks  of  the 
type  under  discussion  is 

rr  u dN' 

p«  ■ “P  |-//dA  ",i  TErj  <”> 

A rather  surprising  relation  exists  between  the  fracture  sta- 
tistics for  volume  distributed  cracks  and  the  fracture  statistics 
for  surface  distributed  cracks.  Let  us  assume  that  uniaxial  data 
are  used  to  determine  the  fracture  statistics  for  specimens  in 
uniform  biaxial  tension.  It  has  been  shown  [36]  that  the  result 
is  the  same  for  both  theories,  i.e.,  it  does  not  matter  whether  we 
use  Eq . (8)  or  Eq.  (37).  This  is  quite  a convenience  for  those  who 
like  to  do  their  own  calculating.  Finding  the  biaxial  stress  sta- 
tistics for  Weibull's  3-parameter  form  using  Eq.  (8)  generally 
requires  the  use  of  a large  digital  computer.  But  using  Eq.  (37) 
the  solution  can  readily  be  found  using  a simple  programmable  hand 
computer  such  as  the  HP-25  or  SR-56. 

There  are  materials  having  surface  distributed  cracks  but  no 
interior  cracks,  e.g.  glass,  but  probably  no  materials  with  volume 
distributed  cracks,  but  no  surface  cracks.  It  is  known  from  frac- 
ture mechanics  that  a crack  near  the  surface  is  generally  weaker 
than  an  identical  crack  far  from  the  surface.  This  complicates  the 
interpretation  of  laboratory  data  on  small  specimens  and  the  appli- 
cation of  the  results  to  larger  specimens.  It  becomes  necessary  to 
separate  out  the  surface  and  volume  contributions  to  the  fracture 
statistics.  An  analysis  of  this  type  has  recently  been  carried 
out  by  Rebholz  and  Teter  [37]. 


MICROSTRUCTURAL  CONSIDERATIONS 


So  far  our  discussion  of  fracture  statistics  has  been  based 
on  a combination  of  extreme  value  theory  and  fracture  mechanics. 

The  crack  density  function  had  to  be  inferred  from  test  data,  and 
it  was  pointed  out  that  extrapolation  of  this  function  to  stresses 
outside  the  region  bounded  by  the  highest  and  lowest  observed  frac- 
ture stresses  is  unwarranted.  Thus  without  additional  sources  of 
information  it  is  not  possible  to  obtain  reliable  estimates  of  the 
fracture  behavior  of  very  much  larger  specimens  or  the  stresses  at 
which  extremely  low  probabilities  of  failure  can  be  assured. 


To  resolve  this  difficulty,  McClintock  [23]  recently  proposed 
a theory  for  the  crack  density  function  based  on  microstruc tural 
considerations.  He  assumed  that  cracks  are  random  aggregations  of 
imperfectly  bonded  pairs  of  grains.  If  the  probability  that  two 
adjacent  grains  are  unbonded  is  q,  the  probability  that  two  such 


unbonded  pairs  are  side  by  side  is  q , etc.  In  this  manner  the 
statistics  of  crack  size  can  be  found,  and  using  the  value  of  the 
critical  stress  intensity  factor  for  the  material  this  can  be  con- 
verted into  statistics  of  critical  stress.  From  this  in  turn  one 
can  obtain  the  statistics  of  failure.  The  material  model  employed 
by  McClintock  was  a 2-dimensional  one  in  which  all  grains  were 
identical  rectangles  and  all  cracks  were  normal  to  the  applied  ten- 
sile stress.  It  was  shown  that  the  probability  of  a crack  being 
longer  than  a is  exp  (-  a/X],  and  that  the.  resulting  fracture  sta- 
tistics do  not  approach  an  asymptotic  form  for  large  specimens. 

This  theory  was  similar  in  many  ways  to  one  published  earlier 
by  Fisher  and  Holloman  [38],  These  authors  considered  randomly 
oriented,  penny-shaped  cracks  in  a 3-dimensional  elastic  solid. 

They  made  the  judgmentally-based  assumption  that  the  probability 
of  a crack  having  a radius  larger  than  r is  exp  [-  r/X].  Using  this 
size  distribution  they  used  fracture  mechanics  to  obtain  a critical 
stress  distribution.  Employing  the  assumption  oe  = On»  they  then 
worked  out  the  fracture  statistics.  The  relation  of  their  work  to 
Weibull's  and  other  WLT's  escaped  general  attention  because  they 
retained  the  binomial  form  in  their  calculations  rather  than  con- 
verting to  exponentials  [see  Eq.  (1)]. 

Batdorf  [39]  applied  the  McClintock  concept  of  the  origin  of 
cracks  to  a consideration  of  randomly  oriented  penny-shaped  cracks 
in  a 3-dimensional  polycrystal.  He  found  that  the  probability  of 
a crack  having  a radius  greater  than  r is  exp  [-  r^/X^].  The  frac- 
ture criterion  employed  was  that  expressed  in  Eq.  (31). 


Among  the  principal  findings  of  this  theory  are  the  following: 

(1)  The  total  number  of  cracks  is  finite  rather  than  infinite 
as  in  Weibull  theory.  The  most  probable  critical  stress  is  sc  - 1 
where  s is  a reduced  stress  to  be  discussed  later.  Weibull's  2- 
and  3-parameter  forms  imply  that  the  number  of  cracks  with  a given 
fracture  stress  increases  monotonically  with  stress;  thus  the  most 
probable  fracture  stress  is  infinity. 


(2)  The  probability  of  failure  in  simple  tension  is  given  by 

(38) 


Pf(s)  = 1 - exp 


- v / (x  - -A  7 exp  (-  1/sc) d sc 
° Sc 


where  N0V  is  the  total  number  of  cracks.  For  a sufficiently  small 
probability  of  failure 
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Pf(s)  ^ (NVq/4)s4  exp  (-  1/s4)  (39) 

(3)  Examining  Eq.  (39)  we  see  that  Pf(s)  = 0 for  s = 0,  and 
Pf(s)  > 0 when  s > 0.  However,  Pf(s)  cannot  be  expressed  as  a 
power  series  in  s because  at  s = 0,  Pf  and  all  its  derivatives 
vanish.  It  turns  out  that  over  a finite  range  of  Pf  such  as 
0.01  5 Pf  < 0.99,  the  prediction  of  this  theory  can  be  fitted 
very  closely  using  Weibull's  3-parameter  form.  However,  m and  Ou 
are  both  very  weak  functions  of  the  volume  instead  of  being  volume 
independent  as  in  Weibull  theory.  There  is  some  experimental  con- 
firmation of  the  predicted  direction  of  change  of  Weibull  param- 
eters with  volume  change  [39]. 

(4)  The  reduced  stress  is  related  to  the  actual  stress 
through  the  equation 

7 _ .0.25/  ,3 

s = 2 0 A / KT  7T 
o / Ic 


£.n  (q  *) 


0.25 


(40) 


where 

A is  grain  cross  sectional  area 
o 

Kjc  is  the  critical  stress  infinity  factor  for  mode  I 

q is  the  probability  that  adjacent  grains  will  be  unbonded. 

Knowledge  of  how  the  fracture  statistics  depend  on  grain  size, 
effectiveness  of  intergranular  bonding,  and  critical  stress  inten- 
sity factor  should  be  very  helpful  in  determining  what  changes  in 
processing  variables  will  effect  desired  improvements  in  material 
properties.  Results  qualitatively  similar  to  those  just  discussed 
are  obtained  with  McClintock's  theory  and  that  of  Fisher  and 
Holloman,  except  that  the  latter  considered  only  a structureless 
elastic  solid,  and  therefore  contained  no  analogue  to  Eq.  (40). 

More  recently  McClintock  has  refined  his  2-dimensional  model 
by  using  regular  hexagons  for  grains  and  assuming  grain  boundaries 
have  a Weibull-type  strength  distribution  [40].  Such  a model 
assumes  no  cracks  are  present  initially,  but  as  the  load  increases, 
cracks  are  created,  grow,  coalesce,  and  eventually  cause  the  load 
capacity  of  the  specimen  to  pass  its  peak  and  decline. 

Another  theory  in  which  cracks  are  created  by  the  loading  pro- 
cess has  been  devised  for  polygraphites  [41  ].  This  theory  is  based 
on  a consideration  of  the  low  tensile  strength  of  a graphite  crys- 
tal in  the  c-direction.  The  grains  are  assumed  to  be  randomly 
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oriented,  so  that  here  and  there  a chance  aggregation  of  grains  of 
nearly  the  same  orientation  are  arrayed  in  the  shape  of  a penny. 

Such  an  array  opens  up  and  acts  like  a crack  when  the  applied 
stress  exceeds  the  capacity  of  all  of  its  constituent  grains. 

Two-dimensional  models  are  generally  of  qualitative  rather 
than  quantitative  value,  but  they  help  guide  our  thinking.  Three- 
dimensional  models  are  obviously  preferable  in  principle,  but 
their  complexity  requires  greater  development  time  and  effort.  At 
present,  simplifying  assumptions  are  used  whose  influence  on  the 
accuracy  of  predictions  is  somewhat  uncertain.  However,  statistical 
theories  of  fracture  incorporating  extreme  value  theory,  fracture 
mechanics  and  microstructural  considerations  are  as  yet  in  their 
infancy.  The  present  writer  believes  that  such  theories  offer  the 
greatest  long  range  promise  for  future  progress. 

STATIC  FATIGUE 

Up  to  this  point  we  have  discussed  primarily  idealized  short 
term  brittle  fracture  in  which  preexisting  cracks  are  unaffected 
by  increasing  stress  until  their  strength  is  exceeded,  at  which 
time  they  expand  suddenly  and  fracture  the  specimen  or  structural 
part.  Some  attention  was  devoted  to  situations  in  which  cracks 
are  created  and  grow  as  the  stress  increases.  We  now  consider 
fracture  resulting  from  subcritical  crack  growth  — slow  growth  of 
cracks  with  passage  of  time,  even  at  constant  stress.  Whereas  the 
objective  in  short-term  fracture  is  to  predict  the  probability  of 
failure  in  a given  stress  state,  in  static  fatigue  it  is  to  predict 
the  time  to  reach  a given  probability  of  failure  when  the  stress 
state  is  specified. 

It  has  been  determined  by  a number  of  investigators  conduct- 
ing mode  I crack  propagation  studies  that  for  a given  system 
[material,  temperature,  environment,  etc.]  there  is  a unique  rela- 
tionship between  crack  velocity  and  the  crack  tip  intensity  factor 
Kj  [42,43].  This  is  usually  expressed  in  the  form 

v *=  A K’J  (41) 


If  a is  the  crack  length 


where  Y is  a geometrical  factor.  Combining  these  relations  with 
Eq.  (41)  it  is  readily  shown  [44]  that  the  time  to  failure  is 
given  by 

T = 2 [Kj~n  - K^“n]  [(n  - 2)  o2  AY2]  (44) 

where  K.  is  the  initial  value  of  KT . 

Ii  I 

The  problem  now  is  how  to  go  from  a laboratory  situation  with 
an  artificial  crack  of  prescribed  size  growing  in  mode  I to  the 
general  situation  in  service,  in  which  an  unknown  distribution  of 
randomly  oriented  natural  cracks  undergoes  mixed  mode  growth  in  an 
arbitrary  stress  state.  Such  a goal  has  yet  to  be  achieved.  In 
fact  it  is  beyond  the  state  of  the  art  in  fracture  mechanics,  since 
it  involves  a knowledge  of  the  growth  rate  and  eventual  critical 
stress  of  a non-planar  crack.  However,  significant  progress  is 
being  made. 

Most  theoretical  work  in  this  area  depends  on  three  basic 
assumptions  or  limitations  (1)  the  stress  state  is  uniform  simple 
tension,  (2)  crack  planes  are  assumed  to  be  normal  to  the  applied 
'stress  both  for  calculating  crack  size  and  for  determining  crack 
velocity,  (3)  the  short-term  fracture  statistics  are  adequately 
described  by  Weibull's  2-parameter  form. 

Davidge,  McLarin,  and  Tappin  [44]  have  used  these  assumptions 
to  develop  strength-probability-time  (SPT)  relations.  An  SPT  dia- 
gram for  alumina  is  shown  in  Fig.  5.  Such  a diagram  can  be  used  to 
find  the  stress  corresponding  to  an  acceptable  probability  of  fail- 
ure during  the  design  life  of  a structural  element. 

Evans  and  Wiederhorn  [45]  have  shown  how  the  statistics  of 
failure  are  affected  by  prior  proof  testing.  The  minimum  time  to 
failure  is  found  by  noting  that  no  crack  longer  than  the  critical 
length  for  the  proof  test  can  be  present  in  any  of  the  surviving 
specimens.  They  then  use  Weibull  statistics  and  Eq . (43)  to 
obtain  the  actual  distribution  of  crack  size  in  the  surviving  speci- 
mens. From  this  they  solve  for  the  probability  of  fracture  as  a 
function  of  stress  and  time  under  load.  Results  for  soda-lime 
glass  in  water  are  shown  in  Fig.  6.  The  comparison  between  theory 
and  experiment  suggests  that  the  former  may  be  conservative  by 
between  a half  and  a whole  order  of  magnitude.  Discrepancies  of 
this  sort  may  be  due  in  part  to  errors  in  determining  material 
parameters  and  propagation  constants.  The  manner  in  which  these 
affect  the  statistics  of  fracture  in  static  fatigue  have  recently 
been  analyzed  in  some  detail  [46,47]. 
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DYNAMIC  FRACTURE 

In  dynamic  fracture  studies,  an  intense  stress  pulse  of  very 
short  duration  passes  through  the  material.  Many  cracks  grow. 

After  passage  of  the  stress  wave,  the  specimen  may  still  be  in  one 
piece  but  in  damaged  condition,  or  it  may  be  fragmented  into  a few 
or  many  pieces. 

In  treating  internal  crack  damage,  Curran,  Shockey  and  Seaman 
i47,48]  have  developed  a fracture  model  that  takes  into  account 
three  aspects  of  the  fracture  process  not  usually  included  in  treat- 
ments of  short-term  fracture:  (a)  nucleation  of  cracks  as  a func- 

tion of  stress;  (b)  growth  of  cracks  as  a function  of  stress,  time, 
fracture  toughness,  and  initial  crack  size;  (c)  decrease  in 
strength  and  stress  attenuation  with  increasing  damage. 

The  nucleation  rate  is  assumed  to  have  the  form 

N - No  exp  [(o  - oj/oJ  (45) 

where  anQ  is  the  threshold  stress  for  nucleation,  and  Nq  and  are 
constants.  The  cracks  are  assumed  to  be  nucleated  with  a distribu- 
tion of  radii  given  by 


Fig.  6.  A Comparison  of  Failure  Times  After  Proof  Testing  [45] 

AN(r)  = ANq  exp  r/rj  (46) 

3 

where  AN0  is  the  total  number/cm  nucleated  in  the  time  interval  t 
to  t + At,  AN  is  the  number  having  a radius  greater  than  r,  and 
is  a constant.  The  growth  of  cracks  was  tested  in  many  ways, 
the  most  successful  being  the  assumption  that 


r = r ( o - cr  ) /4  n (47) 

go 

where  ri  is  the  crack  tip  viscosity.  The  threshold  stress  for  crack 
growth  is  assumed  to  be  given  by 


o *=  0.5  K,  J v It 
go  1c  v 


where  is  the  plane  strain  fracture  toughness. 


(48) 
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As  the  stress  wave  progresses  through  a material,  energy  is 
extracted  in  creating  and  enlarging  cracks,  so  the  pulse  is  atten- 
uated. The  calculation  of  the  number  and  size  of  cracks  as  a 
function  of  time  and  location,  their  effects  on  the  stiffness  and 
damping  of  the  material,  and  the  effects  of  these  in  turn  on  the 
stress  pulse  propagation  is  a strongly  coupled  and  involved  process 
requiring  a computer  code.  Accordingly  we  must  expect  to  be  limi- 
ted to  numerical  results  for  the  particular  cases  chosen  for 
investigation.  However  the  number  and  size  distributions  of 
cracks  in  damaged  specimens  and  fragment  size  in  fragmented  speci- 
mens are  in  reasonable  agreement  with  experiment. 

In  the  case  of  extremely  short  stress  pulses  the  diffraction 
of  the  wave  by  cracks  and  delay  times  in  the  initiation  of  crack 
growth  become  important.  These  and  other  refinements  in  the  theory 
have  been  discussed  by  Kalphoff  and  Seaman  [50]  and  by  Vardar  and 
Finnie  [51]. 


CONCLUDING  DISCUSSION 

It  appears  from  the  foregoing  that  Weibull's  theory  for  uni- 
axial stress  states  is  essentially  correct  as  it  stands,  except 
for  the  limitations  that  it  applies  only  to  tensile  type  fractures, 
and  that  caution  must  be  exercised  in  making  predictions  implying  a 
knowledge  of  n(o)  outside  of  the  stress  range  in  which  it  has  been 
established  by  experiment.  We  have  shown  that  Weibull's  treatment 
of  polyaxial  stress  states  implies  the  assumption  that  cracks  are 
shear  insensitive,  i.e.,  cre  - on.  Eqs.  (8)  and  (11)  allow  for  the 
use  of  arbitrarily  chosen  fracture  criteria  In  analyzing  polyaxial 
stress  states.  Present  evidence  suggests  that  Eq.(31)  is  somewhere 
near  right,  and  leads  to  better  results  than  either  Weibull  theory 
or  the  assumption  of  independence  of  principal  stresses.  More  re- 
search is  needed,  however,  in  the  areas  of  crack  interaction  and 
shear-type  fractures,  where  weakest  link  theory  does  not  apply. 


Theories  including  due  consideration  of  extreme  value  theory, 
fracture  mechanics  and  material  microstructure  i -ve  only  recently 
been  introduced  into  the  literature.  Although  complex,  they  offer 
the  greatest  long  range  promise  and  should  be  developed  much  more 
completely. 

Much  progress  has  recently  been  made  in  the  statistical  treat- 
ment of  static  fatigue,  including  the  effects  of  proof  testing. 

The  accuracy  of  predictions  is  impaired  somewhat,  however,  by  the 
tacit  assumption  that  the  critical  crack  is  normal  to  the  applied 
stress.  Also  little  has  been  done  to  analyze  time  to  failure 
under  polyaxial  stress  conditions. 


In  dynamic  fracture,  considerable  success  has  been  achieved 
in  accounting  for  experimental  data  on  crack  damage  and  fragmenta- 
tion in  the  case  of  materials  that  have  been  studied  fairly  thor- 
oughly. There  are  so  many  material  parameters,  however,  that  opti- 
mizing their  values  to  match  theory  to  experiments  involving  all 
of  them  simultaneously  may  not  result  in  reliable  values  for  each 
parameter.  This  in  turn  makes  transfer  of  the  knowledge  gained  to 
untested  materials  difficult. 

We  conclude  from  all  this  that  a lot  of  progress  has  been 
made  in  statistical  theory  of  fracture,  and  that  useful  techniques 
exist  to  guide  designers  in  their  consideration  of  short-term, 
long  duration,  and  dynamic  fracture  conditions.  It  is  also  evident 
that  a lot  of  work  remains  to  be  done. 
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